NAS A Technical Memorandum 106379 
ICOMP-93-41; CMOTT-93-15 



Time Dependent Turbulence Modeling and 
Analytical Theories of Turbulence 


R. Rubinstein 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
Lewis Research Center 
Cleveland, Ohio 


(NASA-TM-10*'J79) TIME DEPENDENT 
T ) f ... > J ^ f } r C M ) f) | [ M G AND ANALYTICAL 
THi'.i'lh ICS Hr TOR ’ULLNCE (NASA) 

2 n o 


N94-1 7385 


Unc 1 as 


G3/3A 0193175 


November 1993 







TIME DEPENDENT TURBULENCE MODELING AND 
ANALYTICAL THEORIES OF TURBULENCE 

R. Rubinstein 

Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


Abstract 

By sim plifying the direct interaction approximation (DIA) for turbulent shear flow, 
time dependent formulas are derived for the Reynolds stresses which can be included in two 
equation models. The Green’s function is treated phenomenologically, however following 
Smith and Yakhot (Theor. Comput. Fluid Dyn., 4, 197 (1993)), we insist on the short 
and long time limi ts required by DIA. For small strain rates, perturbative evaluation 
of the correlation function yields a time dependent theory which includes normal stress 
effects in simple shear flows. From this standpoint, the phenomenological Launder- Reece- 
Rodi model is obtained by replacing the Green’s function by its long time limit. Eddy 
damping corrections to short time behavior initiate too quickly in this model; in contrast, 
the present theory exhibits strong suppression of eddy damping at short times. A time 
dependent theory for large strain rates is proposed in which large scales are governed 
by rapid distortion theory while small scales are governed by Kolmogorov inertial range 
dynamics. At short times and large strain rates, the theory closely matches rapid distortion 
theory, but at long times it relaxes to an eddy damping model. 
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I. Introduction 


D. C. Leslie 1 proposed solving the direct interaction approximation (DIA) equations 
for shear turbulence 2 by treating the shear terms as a weak perturbation of an isotropic 
turbulent background state. At lowest order, the corresponding perturbation series yields a 
linear relation between Reynolds stress and strain rate in which the eddy viscosity depends 
only on the correlation function and response function of the DIA theory of isotropic 
turbulence. Because DIA is a time dependent theory, this approach leads naturally to a 
time dependent theory of shear turbulence. The goal of this paper is to develop these 
time dependent theories explicitly, incorporating some recent observations of Smith and 
Yakhot 3 on the short time and long time behavior of turbulence. In contrast, Leslie’s 
main concern was the long time, steady state limit of the theory. Evaluating the second 
order terms in Leslie’s expansion leads to a time dependent generalization of Yoshizawa’s 4 
nonlinear eddy viscosity representation of turbulence in which the Reynolds stresses are 
quadratic functions of the mean velocity gradient . 5 Finally, the restriction of Leslie’s theory 
to weak shear is removed by s ummi ng the perturbation series to all orders. For simple 
shear flow in which dUi/dxj = SSuSj?, the summation is accomplished with the help of 
rapid distortion theory (RDT ). 6 

The present theory of weakly sheared turbulence can be compared to a standard 
phenomenological time dependent model, the Launder- Reece- Rodi (LRR) Reynolds stress 
transport model . 7 From the viewpoint of the present theory, the LRR model arises by 
incorrectly replacing the DIA Green’s function by its long time limi t.. This simplification 
misrepresents the short time response of turbulence to shear. An important consequence is 
that eddy damping corrections to short time RDT behavior initiate too quickly in the LRR 
model. This causes excessive growth of turbulence kinetic energy at short times in highly 
strained homogeneous shear flow. The present theory exhibits a strong suppression of 
eddy damping at short times. It therefore also predicts a much wider frequency range over 
which RDT correctly describes oscillating shear flow: in the present theory, the corrections 
to RDT are of order u >~ 2 where u> is the oscillation frequency, but are of order u > -1 for the 
LRR model. This may explain the success of RDT based models in computing oscillating 
flows . 8 

The present proposal for strongly sheared turbulence leads to a ’’two scale” picture of 
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shear turbulence in which large scales are rapidly distorted while small scales axe governed 
by Kolmogorov inertial range dynamics. The theory can be described either as RDT with 
generalized eddy damping, or as RDT with a modified total strain.The introduction of a 
phenomenological modified total strain to improve the agreement of RDT with experiments 
has been advocated in the RDT literature ; 6 ’ 9 in the present theory, the modified total strain 
is deter min ed by the Green’s function of isotropic turbulence. This theory is intrinsically 
time dependent; however, it can be simplified by assuming that turbulent states at any 
strain rate can be long time limits. Then for simple shear S = dJJ\ j dx 2 held for infinite 
time, the theory formally reduces to a relation 
K 2 

T12 = C v (r))—S 

with 10 77 = SKje. There are analogous formulas for the normal stresses. Phenomenological 
expressions of this type have been proposed , 11 but in the present theory, C v (r)) is exactly 
determined by RDT. 

II. Formulation of the Theory 

A. Simplified DIA Analysis of Shear Turbulence 

Leslie’s theory of shear turbulence 1 can be derived from the generalized Langevin 
model for isotropic turbulence 12 

d ft 

— + vk 2 ) «i(k,t) -f / ds rj(k,t,s) Ui(k,s) = /i(k,t) 
ot 

in which Ui is the random velocity field, r/(fc,t, s) is a deterministic eddy damping factor, 
and / is a random force. A Fourier space representation is used and k 2 = k • k. The 
properties of r)(k,t,s) and / axe given in detail in Ref. 12. They depend on the two 
time correlation function of the velocity field, so that the linearity of this equation is only 
apparent. Eq. (1) is also a generic model in the statistical mechanical theory of transport 
coefficients . 13 DIA gives the exact correlation function for a suitable model of this type . 12 

Suppose that some external agency, such as shear or buoyancy forces, is present. We 
will generalize Eq. (1), representing the effect of the external agency by adding a suitable 
force Fi to the right side: 

+ vk 2 ) ui(k,t) + J ds r?(M,s) «i(k,s) = /i(k,t) + -F;(u p (k,f)) ( 2 ) 
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This general model, which can perhaps be attributed to Leslie, is a simplification for 
each force F of a corresponding complete DIA theory. It simplifies the dynamics by 
ignoring any effects of the external agency on either the eddy damping or the random 
force; in particular, the eddy damping always remains isotropic. Except for some tentative 
proposals of Tchen 14 and the suggestions of Cambon et al, 15 such effects have been little 
investigated. Ignoring them amounts to treating shear turbulence, for example, as the 
outcome of straining by large scales acting against isotropic eddy damping. Although 
this picture is oversimplified, it is a plausible starting point and should capture some of 
the physics of shear flow. Modifying the damping in the region of strong effects of F may 
overcome some of the limitations of the model. 16 In any case, a more complete investigation 
based on a full DIA analysis will be considerably more difficult. 

With these assumptions, the model equation for shear turbulence is 

ds T](k,t,s) Uj(k, s) = + S ip (k,t)u p (k,t) (3) 

where 

Sim = Aim 2fe k{k p A prn -t - 8i m k)A ir d j dk r 
Aim = dUi/dx m 

Note that the viscous term of Eq. (2) has been ignored in Eq. (3). The shear terms 
result 17 from Reynolds averaging the Navier Stokes Equations with a mean velocity field 
Ui = AijXj. Reynolds averaging also introduces the convective derivative in Eq. (3) which 
replaces the time derivative in Eq. (2). Since convection by the mean flow does not affect 
eddy damping or straining, this effect will be ignored in what follows. Closely related 
models have been proposed and investigated by Cambon et al 15 in the context of EDQNM 
especially for rotational effects. Since our goal is to derive single point models rather than 
to study spectral dynamics, Leslie’s much simpler formulation seems adequate. 

Eq. (3) can be rewritten in terms of a Green’s tensor or response tensor 

Gij( k,t,s) = G(k,t,s)Pij(k) 



where 

-Pij(k) = &ij ~ k{kjk~ 2 
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as 


( 4 ) 


( 5 ) 


( 6 ) 

In what follows, the condition Eq. (6) on the Green’s function will always be understood. 
Moreover, in order to avoid additional notation, the symbol G will be used to denote 
different Green’s functions dependent on different numbers of arguments. Eq. (3) will be 
applied to derive Reynolds stress models, however it should be possible to deduce K and 
e transport equations from it as well. 

B. Relation to RDT 

Eq. (3) with both the random force and the eddy damping terms neglected is just 
RDT. 6 An approximate condition that RDT apply is therefore that 

■q(k,t,s)ds « ( ApqApq y! 2 (7) 

Because 77 models eddy damping due to nonlinear interaction, this condition agrees with 
the usual idea that RDT applies when the shear dominates nonlinearity. At short times, 
the left side of Eq. (7) is 0(t)\ therefore, the short time response is always governed by 
RDT in this theory. At finite times, Eq. (7) states a scale dependent criterion for the 
applicability of RDT which will be discussed later. 

In the RDT limit with Tj(k,t,s) = 0, Eqs. (5), (6) become simply G(k,t,s ) = 1. Thus, 
a second criterion for the applicability of RDT is 

|G(fc,i, s) — 1| << 1 (8) 



Ui(k,t) = f ds G(k,t,s) [/j(k,s) + P ig (k)S, p (k,s)u,(k,s)] 

Jo 


where G satisfies the integrodifferential equation 


dG(k,t, s)/dt + J dr r](k,t,r) G(k,r,s) = 0 


and the conditions 


G(k,s,s) = 1 
G(k,t,s) = 0 for t < s 
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In what follows, eddy damping will be described by G rather than by 77 , and Eq. (4) will 
be used as the basis of the theory. 

C. Stationary Green’s Functions 
The time stationary case 


G(k,t,s) = G(k,t-s) (9) 

should be characterized by universal inertial range forms with Kolmogorov similarity 

G(k, t-s) = G(0, C = e i/ 3 jfc 2 / 3 ( t _ ( 10 ) 

No attempt to calculate the time dependent function G(k,t — s ) theoretically has yet 
succeeded. DIA itself is inconsistent 18 with the Kolmogorov scaling of Eq. (10); the 
Lagrangian modification 19 of DIA which restores Kolmogorov scaling does not give sat- 
isfactory predictions for the long time behavior of time correlations . 16 At this time, it is 
therefore necessary to postulate functional forms for <?(£). However, Smith and Yakhot 
observed 3 that the short and long time limits of G alone have important consequences. It 
follows from Eq. (5) that at short time separations, 

G(k,t — s) = 1 + 0((t — s) 2 ) for (t — s) ~ 0 ( 11 ) 

and it is generally believed that at long time separations, the function (?(£) decays expo- 
nentially, 

G(k,t — s) ~ exp(— Co£) for (t — s) ~ 00 ( 12 ) 

with Cd a universal constant: in the Yakhot-Orszag theory , 20 Cd — -49. This limi t, 
corresponds to eddy damping. We will follow Smith and Yakhot 3 by leaving the functional 
form of Cr(f ) unspecified, but insisting on the limits in Eqs. (11), (12). 

An important observation is that the long time limit, Eq. (12), does not satisfy the 
short time limit Eq. (11). In fact, reference to Eq. (5) shows that the long time limit is 
defined by the singular damping function 

T](k,t,s) = — s) (13) 
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which. ’’Maxkovianizes” Eq. (3) as in Kraichnan’s test field model. 12 It will be shown later 
that the LRR stress transport model also assumes Markovian eddy damping. 

A consequence of the short and long time limits follows from the definition of £ in Eq. 
(10): at any fixed nonzero (t — 5 ), £ is small for large scales and large for small scales. Eqs. 
(8), (10) imply that sufficiently large scales are governed by the short time limit, RDT, while 
sufficiently small scales are governed by the long time limit , inertial range eddy damping . 
This observation suggests a ”two-scale” theory of shear turbulence. 

D. Nonstationary Green’s Functions 

For completely general conditions, the time dependence of G(k y t , s) can only be found 
from DIA. In the present simplified theory, the time dependence must be postulated in- 
stead, recognizing that some form of universality is indispensable in turbulence modeling. 
This will restrict the applicability of the theory, but such restrictions are inevitable in 
any case: there are time dependent problems accessible to DIA, such as the relaxation of 
turbulence with strong k space anisotropy, or the generation of a Kolmogorov spectrum 
from an arbitrary initial spectrum, which cannot be usefully described at the single point 
level. 

Let the inertial range be characterized by its time dependent dissipation rate e(t) and 
inverse integral scale fco(^)- In the Yakhot-Orszag theory, 20 k$ is defined so that 

K(t) = f E{k,i)dk 

J fco(t) 

where K denotes the turbulence kinetic energy; for the Kolmogorov spectrum written as 

E = for k > ko (14) 

where Ck is the Kolmogorov constant, K and &o are related by 

K(t) = |c KE (t) 2/3 t 0 (()- 2/3 (15) 

Either pair of functions e(t) and k 0 (t) or e(i) and K(i) defines a time dependent inertial 
range: the second pair will be assumed to be known from the solution of a two equation 
model. 
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Define the frequencies 


6(k,i) = e 3/3 (t)k 2/3 for k > k a (t) 

m = e(i )/m 


A universal time dependent damping function can be postulated by ass uming that a sta- 
tionary Green’s function (?(£) satisfying the limits Eqs. (11), (12) is known, and replacing 
the similarity variable £ = 9(t — s ) of Eq. (10) by the generalization J* 0 dr, so that 


G{k,t, 



9{k,r)dr) 


(16) 


This postulate has two consistency properties: it reduces to the stationary form when the 
inertial range is time independent, and it is exact for the singular case Eq. (13) in which 
G(0 = e~ c °t. 


A simpler formulation, closer in spirit to single point modeling, is to treat damping 
as scale independent by setting G(k,t,s ) = G(t,s) only. In this case, the damping is due 
to the action of the inertial range as a whole. It will be convenient to call this type of 
model a global damping model. The appropriate similarity variable in the stationary case 
is (t — s)0. The short time limit is 


G(t- s) = 1 + 0((t - s) 2 ) for (t - s) ~ 0 (17) 

and the long time limit is 

G(t - s) ~ e -c R ('/K)(t-s) for (< _ ^ oo ( 18 ) 

where Cr is another universal constant: in the Yakhot-Orszag theory, 20 Cr ~ 1.58. As- 
suming that G{t — s ) satisfying Eqs. (17), (18) is known, we can define by analogy to Eq. 
( 16 ). 

G(t,s) = G(f Q(r)dr) (19) 

Dropping the k dependence in Eq. (5) shows, by analogy to Eq. (13) that the long time 
limit corresponds to the singular damping function 


= ®(s)8(i — s ) 


( 20 ) 
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Eq. (19) is exact in this case, for which G(t - s) = exp(— Cji(i — «)©)• 

III. Time Dependent Eddy Viscosity 

To derive the t im e dependent eddy viscosity, we follow Leslie 1 and expand Eq. (4) in 
powers of the mean strain rate about an isotropic background state u^: 

u = u (°) +U W + ... ( 21 ) 

We have assumed that the force / is independent of the mean strain rate; therefore, the 
effect of / is absorbed entirely in the background state, and is given by 

,i) = t ds G(k,t,s) Pim{ k) S mn (k ,s) u< n °> (k, s) (22) 

Jo 

If, corresponding to Eq. (21) the single time correlation function is expanded as 

Q = Q ^ + - • * (23) 

where 

<?Sf(M) = < (k,<) uf\-k,t) + «i° ) (k,t)rij 1) (— k,t) >/6(k) 

and the higher order correlations are defined similarly, then Eq. (22) implies 

<?im(M) = / ds G(k,t,s)(—Ai r + 2kik p k- 2 A pr )P rm Q (0 \k,t,s) 

Jo 

+(*m) + G(k,t,s)* r A rn Qj- Q { £{k,t, S ) (24) 

— s)k r A rn G(k,t,s) 

where ( im ) denotes index interchange in the immediately preceding term. In view of 
the isotropy of the background field, Q^(k, f,s) = Q(°)(fc,t,s)Pij(k). The occurence of 
two time correlation functions in the formula for the single time correlation function is 
characteristic of DIA. The time stationary form of this equation was stated by Leslie. 1 We 
follow Smith and Yakhot 3 and assume the nonstationary fluctuation dissipation relation 

Q (0) (k,t,s) = Q(°\k,s)[G(k,t,s) + G(fc,M)] 
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A decomposition of the Reynolds stress follows from Eq. (23): 


— (0) + t£> + 


T ij = T .j 


where 

rjfw = I dk «S;>( k,(> 

Note that the sign convention established here is 

Tij = + < Ui Uj > 


Evaluation of the angular integrals in Eqs. (24), (25) leads to 

T ij\i) = ~ J Q ds J Q dk G(k,t,s) 2 

~~ G(k,t,s)k ^ G(fc,t,s)} E(*,s) 5y(*) 
where E(k,s ) is the energy spectrum at time s and 


Sij 


Ay + Ajj 


( 25 ) 


(26) 


We adopt the viewpoint of the Yakhot-Orszag theory 20 and evaluate Eq. (26) over 
the inertial range k > ko only using the similarity form Eq. (16) for the Green’s function 
and the Kolmogorov spectrum, Eq. (14). The result will have the general form 

T iP = / ds T(t,s)Sij(s) 

Jo 

in which the integral kernel T itself depends in a complicated manner on the evolution of 
the inertial range parameters K and e for times between 0 and t. In the stationary case, 
the result can be rewritten as 


= -Ck 

*00 

I d.K «-*/* 

(|C r r)»/J 



dr Sij(rK/e) x 
- jg k 2 / 3 G(« 2 / 3 )G’(« 2 / J )} 


(27) 


Recall that G here denotes the inertial range similarity form Eq. (10). Given a functional 
form for G, the second integral in Eq. (27) is a universal function of the dimensionless 
time variable r. 
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As Smith and Yakhot emphasize, 3 the universal short and long time limits for G imply 
universal short and long time limi ts for Namely, substituting the short time limit 

for G given in Eq. (11) in Eq. (26) leads to 


K = -77 KtS, 


'tj 


15 


*7 


(28) 


in agreement with Crow’s 21 RDT calculation. To derive the long time limit, use the 
stationary form Eq. (10) for G(k,t — s ) and define constants C\ and C 2 by 


f 


G(k,a) 2 da — C\/Q 


s: 


G(k,a) k dG(k,a)/dk da = C 2 /0 


Then ass umin g constant strain rate and substituting in Eq. (26), using Eq. (15) to 
elimin ate kg in favor of K, 


where 


r - -C K * S- 


= ^)/45C K 


There are some other useful special cases of Eq. (26). By setting G = 1, we obtain 
RDT expanded to first order in the mean strain rate. The possibility of a ’’viscoelastic” 
representation of turbulence, as in Eqs. (26), (27) was suggested by Crow 21 and others. 
Crow’s theory arises in the present formalism by assuming viscous damping 

G{k,t,s) = e -^ 2(t - a) 


instead of eddy damping. 

Global damping models greatly simplify these formulas while retaining the idea of 
eddy damping. The global damping analog of the general nonstationary model Eq. (26) is 

T§> = f is G(t,sf K(s) Sij(s) (29) 

An important special case of this formula is obtained by substituting the exponential form 
for G of Eqs. (18),(19); although this form does not satisfy the short time constraint Eq. 
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(17), it will connect the present theory with stress transport models. In this case, Eq. (29) 
becomes 



\exp 


t- 


2Cr®{t) dr K(s) Sij(s) 


Equivalently, 

4 f) = - Ti KSii {30) 

or for constant strain in simple shear in which Aij = SSnSj 2 


2 K 2 /„ . „„ e 

Tl2 = "I 

In the Yakhot-Orszag theory 20 , 2/15 Cr = C v ~ .08 exactly equals the usual eddy 

viscosity constant. 


Comparison with the LRR model 7 requires that both convection and diffusion be 
ignored: these are inhomogeneous effects extraneous to the present analysis. Writing T-j 
for the deviatoric, or anisotropic part of the stress, the simplified LRR model is 


dt l K ij 15 

, p r D , 


, t D dUp D ou 

+ CllTi >9^ +r »di 


,dUi 

ip dx p 

au p 


--ft t D 
3®*i T 5J> 

- -ft t D 
3 ®*J T 9P 


(31) 


Solving this equation by a perturbation series in the strain rate analogous to Eq. (21), 


d (i) . (2) . 

Trf = T- ■ 4- T- V - ’ 4- 

tj * 3 ~ ' 


we find that rjp satisfies Eq. (30) with 2 Cr — C\ . This establishes a simple connection 
between the present theory and LRR, namely that to lowest order in this perturbation 
theory, the LRR model is a global damping model defined by the Green’s function 


G(t,s) = exp f — (7i0(r) dr 


(32) 


which does not satisfy the short time limit Eq. (17). Instead, this choice of Green’s 
function corresponds to the long time limit Eq. (18), to the singular eddy damping of Eq. 
(20), and to Markovian damping in Eq. (3). 
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The LRR model satisfies the Crow constraint on short time behavior, Eq. (28). The 
short time expansion G(t — s) ~ 1 + 0(t — s) of the LRR Green’s function Eq. (32) implies 
that corrections to the Crow constraint are of order t 2 ; Eq. (17) implies a correction of 
order t 3 . Summarizing, the short time corrections to the Crow constraint are 

r ~ t + 0(t 2 ) LRR 

r~ t + 0(t 3 ) present (33) 

r~ f + 0(f 5 ) RDT 9 

These distinctions are important in oscillating shear flow. Consider oscillating simple shear 
flow in which S = dU\jdx 2 is the only nonzero mean velocity component, and let r = 
be the shear stress. Suppose that S is oscillatory: 

S = da/ dt = aiu>e lwt 

In rapidly oscillating flow, it is reasonable to ignore oscillations in K and e and to replace 
them by their time averages. Then © = e/K is constant. Damping is then described by 
the stationary Green’s function G((t — s)0). Then r is given by 

r(t) = ——aiKui [ ds G{Q(t - s)) 2 e™' 

15 Jo 

The steady state solution when u> is large satisfies 22 

r ~ — 4aJf/15 A 0 + -A + . 2 2 , H ]e lu>t 

where 

Ao = G(0) = 1 
A j = 2G(0)G'(0)© 

A 2 = [2G(0)G"{0) + 2G'(0) 2 ] 0 2 

indicating as expected that in the limit u — > oo, the response is elastic and is governed 
by RDT. In the LRR model, G'(0) ^ 0; therefore the corrections to RDT are of order 
u> _1 . But if G satisfies Eq. (17), then Ai = 0 and the correction to the RDT solution 
depends on u>~ 2 instead of on a? -1 . Moreover, the phase lag between stress and strain is 
of order u> -3 . This may explain the success of RDT in solving high frequency oscillating 
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flows. 8 It appears that the LRR model will overpredict the effects of eddy damping at high 
frequencies. 

The corrections to RDT at short times summarized by Eq. (33) are also important 
in transient homogeneous shear flow. Eq. (33) indicates that the short time corrections in 
both LRR and the present theory are associated with the initiation of eddy damping. It 
is known 23 that at high strain rates, the LRR model predicts a much too rapid growth of 
K. Energy growth is due to the onset of turbulence production by eddy damping. In view 
of Eq. (33), eddy damping corrections are of order t 2 in LRR, but are of order t 3 in the 
present theory. This strong short time suppression of eddy damping suggests that RDT 
will apply in the present theory for longer times than it does in the LRR model. This may 
improve the agreement with DNS studies 24 of highly strained turbulence. 

IV. Second Order Analysis 

Calculation of Leslie’s expansion to second order is lengthy but routine. The result 
has the form 


eS?(k,i)= E 

1<2V<6 

+ c WA„(s)A ir (r) + <t N >A ri (s)A rj (r) 

+ e(' f >S i ,A„(s)A, r (r) + + (y) 


in which the are integral operators 


I* 1 ) = f ds G(k,t,s) I* 
Jo Jo 

I (2) = [ ds G{k,t,s) [ 
Jo Jo 

i< 3) = [ ds G(k,t,s) r 
Jo Jo 

/ (4) = f ds G(k,t,s) [ 
Jo Jo 

J (5 > = f ds G{k,t,s) f 
Jo Jo 

j(6) = f ds G(k,t,s) [ 
Jo Jo 


dr G(k,s,r)G (k,t,r)Q m (k,r) 
dr G(k,s,r)k-^[G(k,t,r)Qm(k,r)) 
drG(k, s ,r)k 2 ^[G(k,t,v)QW(k,r)} 
dr G(k,t,r)G(k,\ s,r \)Ql°>(k,r) 
dr G(k,i,r)k-^{G(k,\.,r \)QU\k,r)] 
dr G(k,t,r)k* -£-{G(k,\ s,r |)0<“>(*,r)] 


(33) 


( 34 ) 
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a,( N \...f( N ^ axe the following geometric constants: 



(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

105a( N ) 

27 

-1 

-2 

19 

2 

1 

2 

-1 

105 bW 

20 

6 

-2 

6 

3 

-1 

105c< N ) 

-15 

-15 -2 

5 

2 

_ 15 
2 

-1 

105d<"> 

20 

-8 

-2 

-22 

-4 

-1 

105e( N ) 

10 

24 

6 

3 

12 

3 

105/W 

-4 

24 

6 

10 

12 

3 


|) = G(JW) + G(fc,r,s) 

As in all derivations of this type 4 ’ 5,9 , these constants arise from integrating even order 
products k{kj , .... over spheres k = constant. A useful constraint on the calculation is that 
it must reduce to the rapid distortion results of Maxey 9 when G = 1. Single point models 
are derived exactly as in Sect. II and will not be described explicitly. 


In a global damping model, the integrals in Eq. (34) satisfy the relations 

/(2) = —3 1^ 

I (3) = 12 1 (1) 

j( 5 > = -3I (4) 
j(6) _ 12 j(4) 

and Eq. (33) reduces, after perfor min g the wavenumber integrals, to 

T ij\t) = f ds <?(*,«) f dr G(s,r)G(t,r)K(r) x 

Jo Jo 

{ + lS ip (s)A ip (r) + S jr (,)A ip ( r)] (35) 

24 

— [>S'tp( s )Apj(r) + 5jp(s)A p i(r)] 

32 

+ — — 5tj[A pg (s)A 9P (r) + -4.p9( s )-^-pg( r )]} 


To compare this calculation with LRR, evaluate the second order solution of Eq. (31): 

T ij\t) = ! ds G(t,s) [ dr G(s,r)K(r) x 
Jo Jo 

+ C , 2[Sip(s)Ajp(7’) + S’ J p(s)Ajp(r)] 

+ C , 3[S»p(s)Apj(r) + 5j p (s)Apjfy)] 
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where G is the Green’s function for the LRR model, Eq. (32). The occurence of an 
additional Green’s function in Eq. (35), which arises from the two-time correlation, means 
that at second order, the present theory does not reduce to LRR. However, as in the 
discussion of shear stress, the important difference between the theories is the correction 
to short time behavior: in both the LRR model and the global damping model, the normal 
stresses are of order t 2 at short times. In the LRR model, the corrections are of order t 3 , 
whereas in a global damping model in which G satisfies the short time limit Eq. (17), the 
corrections axe of order t 4 . Thus, the present theory resembles RDT longer in transient 
homogeneous shear flow and at lower frequencies in oscillating shear flow than the LRR 
model. 

V. Models Valid for Large Strain Rates 

The perturbative derivation of these models limits their applicability to moderately 
strained flows. A theory valid at arbitrary strains can be derived by summing Leslie’s 
expansion 1 to all orders. The summation is simplest in the time stationary case which will 
be considered first. In an obvious operator notation, Eq. (22) can be written as 

u (1) = GSu(°> 

where G denotes time convolution by the (stationary) Green’s function, 

G(t — s)a(s) ds 

(only the time arguments have been shown explicitly), and S denotes the action of the 
strain dependent terms. In view of Eq. (3), 

u = GSu( n_1 > = (GS) n tt<°> 

where operator products are understood; therefore, 

u = [I + GS + (GS) 2 +...]u<°> (36) 

and the problem is to find a useful representation for the Neumann series on the right side. 
When G = 1, G is simply time integration, which can be written 

(Go)(i) = ( H * a)(t) = f a(s ) ds 

Jo 


(Ga)(t) = (G * a){t) = [ 

Jo 
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where H denotes the usual unit step function. In this case, the sum in Eq. (36) defines 
RDT, for which the sum can be given explicitly for some important special mean velocity 
gradients Aij. 

Consider a global damping model G(k,t , s) = G(t, s), so that in each term of the series 
in Eq. (36), the k-derivative in S does not act on G. RDT can be written as the special 
case of Eq. (4), 

Ui(k,t) = f ds H(t - s ) P iq (k)S qp {k,s)u p (k,s) (37) 

Jo 

Following Ref. 9, write the solution of RDT for simple shear 

Aimi.i') = 
as 

= M ip (m(k,a(i)),a(<))u p (k,0) 

where c*(f) is the total strain 

a(t) = [ S{s) ds = (H * S)(t) 

Jo 

and m is defined in Ref. 9. Define the modified velocity gradient A* m by 

A; m = X * A im (40) 

where componentwise convolution is understood and the function X is chosen so that 

H*X = G (41) 

therefore 

H * A* m = G * A im (42) 

Comparing Eqs. (37), (40), and (42), it is evident that the solution of Eq. (37) with the 
modified velocity gradient A* m is the solution of Eq. (4) with an arbitrary Green’s function 
G. 

The solution of Eq. (41) is 

X = G’ + 6 (43) 


(38) 

(39) 
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where G' is the derivative of the stationary Green’s function G(t — s) with respect to its 
argument. For simple shear, Eq. (40) reduces to 

5* = X * S (44) 

Therefore the solution of Eq. (4) for simple shear is the modification of Eq. (38), 

Ui(k,<) = M ip (m(k,a*(t)),a*(t))up(k,0) (45) 

where, in view of Eqs. (41), (44), a* is defined by 

a* = H*S* = G*S 

Explicitly, 

a*(f) = f G((t - *)©)$(«) ds (46) 

Jo 

Eqs. (45) and (46) solve the problem of summing Leslie’s perturbation theory in the special 
case of simple shear, a* may be called the modified total strain. 

The introduction of a phenomenological modified total strain has often been suggested 
in the RDT literature 6 ’ 8 ’ 9 as a way to improve the agreement between RDT and data from 
flows which axe not evidently rapidly distorted. The short and long time properties of Eq. 
(46) are interesting from this viewpoint. At short times, G ~ 1 and Eq. (46) becomes 

Q * rv (X ~ S(0)t for t 0 (47) 

indicating that the short time limit of this theory is RDT. At long times and constant 
strain rate S', 

a*~ / e~ c * @a Sds = — fort~oo (48) 

Jo Gr e 

In shear flows nearly in a production equals dissipation steady state, SK/e ~ 3.0; for flows 
which evolve to this state, the modified total strain saturates after growing linearly for 
short times. In the RDT solution for constant strain rate, the total strain grows linearly 
for all times. The phenomenological modifications of RDT suppress this growth by forcing 
saturation of the modified total strain. Here, this saturation is a consequence of Eq. (46). 

Some obvious generalizations of this theory are to nonstationaxy problems and to gen- 
eral wavenumber dependent damping models. Arguing by exact analogy to the stationary 
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case, to develop a nonstationary theory, it would be necessary to interpret the operator G 
as a Volterra operator 

(Ga)(f) = f G(t,s)a(s) ds 

Jo 

and to define the modified velocity gradient as the solution of a Volterra integral equation. 
In view of the phenomenological character of the nonstationary Green’s function, it seems 
preferable just to generalize Eq. (46) directly as 

«*(<)= f G(t,s)S($) ds 
Jo 

When the stationary Green’s function is fc- dependent, the summation fails because 
the k derivatives in the strain operator S also act on G(k , t — s). The direct generalization 
of Eq. (46) for simple shear flow, 

a*(k,t) = f G(k,t — s)S(s) ds (49) 

Jo 

therefore does not lead to an exact s umm ation of Eq. (36). But inserting the wavenumber 
dependent modified total strain of Eq. (49) in Eq. (45), 

«i(k,i) = Mi p (m,k,a*(fc,t)),a*(A:,<))up(k,0) (50) 

does at least sum terms of all orders in the strain rate in Eq. (36); it is therefore a plausible 
approximate summation of this series. Eq. (50) brings about a generalization of the short 
and long time limits of Eqs. (47) and (48): since G(k, t — s) ~ 1 for sufficiently large scales, 
whatever the value of (t — s), the RDT solution for large strain rates can apply even at finite 
times to sufficiently large scales. This idea is also stated in a s imil ar context by Cambon et 
al. 15 However, for sufficiently small scales, G rapidly assumes its asymptotic eddy damping 
form. Eq. (50) states a ”two-scale” theory of shear turbulence in which large scales can be 
governed by RDT while small scales exhibit the eddy damping characteristic of a steady 
state Kolmogorov inertial range. 

V. Approximate Theory of Highly Strained Flows 

There have been recent proposals 11 to modify the standard eddy viscosity 

vt = C V K 2 / e, C v ~ .09 (51) 
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to improve its behavior for large strain rates. These proposals all introduce functions 
C v (t}) in Eq. ( 51 ), where in simple shear flow with Aij = S 8 n 6 j 2 , rj is the dimensionless 
strain rate tj = SK/e. This idea can be attributed to Yakhot et al, 10 who replace the 
constant C € i of the two equation model by a function C^(r)). Applied to the nonlinear 
eddy viscosity model of Yoshizawa, 4 this approach suggests 

2 K 2 

Tij = ~K 8 l3 - cm—Sh 

O £ 

K 3 \ 

+ C r rl( 7 ?) ^2 l^ipAjp — ( 52 ) 

K* 2 

C r 2 {jf) [AipApj + AjpApi Q^tjAqpApq] 

3 ^ 

+ ^ [ApiApj — J^fiijAqpAqp] 

This equation, with constant C V ,C T \,C T 2, and C T 3 cannot be applied to flow regions, like 
the near wall, in which 77 is large: it predicts that the stress ratios Tjj fK all increase with 77, 
completely contrary to the data, and even predicts that some of the normal stresses become 
negative. These problems can be overcome if the functions of 77 in Eq. ( 52 ) are chosen 
suitably. An early, and especially interesting model of this form was proposed by Pope 25 
who solved the ’’algebraic” form of the LRR model due to Rodi 26 for two dimensional 
mean flow. 

A model of this type can be derived for simple shear flow under the additional hy- 
pothesis that Eq. ( 48 ) defines a long time li m it for any value of 77. It must be stressed 
that this is an additional assumption; it is not a consequence of this theory. Analytically, 
this assumption states that a* in Eq. ( 45 ) can be replaced by its formal long time limit, 
7 i/Cr from Eq. ( 48 ), so that 


ui(k,t) = Afj p (m(k, 77/ C#), 77/ CR)u p (kj 0) 


Then 


K 2 

r = -CM — S 

E 


where t — t \ 2 and the function (77) is found from RDT as follows. By forming moments 
and integrating over wavenumbers, Maxey 9 presents t/K graphically as a function of 
a = St, say 

t/K = F{a) 
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Replace a by tj/Cr as in Eq. (53); the result is 


r = kf(t]/Cr) = -cm —s 


where 


c. M = 


F(n/C R ) 


( 55 ) 


and the function F is known from RDT. Note 9 that since F — > 0 when 77 — > oo 


vCv(v) V ->• 00 


in agreement with the observation that strong shear suppresses the shear stress so that 
t/K — * 0 when 77 is very large. Strain dependent coefficients C r i(rj),C T 3 (rj) for Eq. (52) 
are also easily obtained from Maxey’s RDT results. 

This theory helps explain a curious feature of homogeneous shear flow data: the ratio 
t/K is about the same both in fully developed homogeneous shear flow in which 77 ~ 4.4 
and in simple shear flows in energy equilibrium in which 77 ~ 3.0. Therefore, if the usual 
formula Eq. (51) is written as 

t/K = —C v t) 

and is calibrated for equilibrium shear flows, it will predict a viscosity which is too large 
in homogeneous shear flow. The data is summarized in Table I and compared with the 
predictions of Eqs. (54)-(55). The theory predicts that t/K has approximately the same 
value in both flows because in both flows 77 j C r is near the maximum, at about a = 2 , of 
F(a) according to Maxey’s RDT calculation . 9 

Table II compares t/K in near wall channel flow with the theoretical predictions. 
Making this comparison invokes the suggestion of Lee et al 24 that near wall turbulent states 
are similar to highly strained homogeneous shear flow. As in the homogeneous shear flow 
comparison, there is qualitative agreement with the trends, but the theory underpredicts 
the reduction in t j K as 77 increases. In both Tables, results for the theory with transport 
corrections are also listed. These arise as follows: transport effects (convection by the mean 
flow and turbulent diffusion) are inhomogeneous effects which have been ignored in this 
theory. Rodi’s 26 algebraic transport correction is useful because it is model independent; a 
rough way to incorporate it in this model is to reduce the viscosity by the factor C\/{C\ + 
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P/s — 1 ) where C 1 is the constant of the LRR model Eq. (31) and P is production. We 
have set C\ = 1.6 in applying this correction. It improves the agreement with the data in 
some cases and greatly overcorrects in others; however, the comparisons axe not me aningf ul 
unless some correction for transport is made. 

A referee has suggested comparison with Pope’s explicit solution 25 of the algebraic 
LRR model. Further comparisons with this model are given in Sect. VI. The comparison 
shows that the algebraic LRR model also reduces the ratio r/K as i 7 increases. The 
constants of Ref. 25 have been used to make the comparisons. The quantitative agreement 
could be improved by adjustments of these constants. 

The quantitative limitations of the theory can be attributed to its basic assumption 
that the large 77 states can be considered long time limits. It is more likely that highly 
strained states are transient; although near wall flow can be steady in time, this steady 
state is maintained by the continual diffusion of highly strained turbulence away from the 
near wall production region into to bulk of the flow. From this (Lagrangian) viewpoint, 
turbulence is highly strained only for a finite time . 16 The simplicity of models like Eq. 
(52) makes them attractive, but they all assume the steady state character of highly 
strained states. Therefore, they should not be accepted uncritically, although their value 
in "regularizing” the stresses in large 77 flow regions will make them useful in calculations. 

In Sects. IV and V, simple shear flow has been emphasized because an explicit RDT 
solution exists for this flow. In principle, the constructions of these sections can be gener- 
alized to any mean velocity gradient, but they will be explicit only when a corresponding 
RDT solution is known. Thus, a theory of this type applicable to near wall calculations in 
square duct flow will require solving RDT for a mean velocity gradient with the structure 


dUJdxi dUj/dx 2 0* 
dU 2 /dx a dU 2 /dx 2 0 
dU z /dx 1 du 3 /dx 2 0 


dUi dU 2 

dxi dx 2 


0 


It should not be assumed that superposing results for simple shear flows will provide a 
good approximation, or that the results will depend on a single parameter like 77. The 
development of a more general theory, in which a rotational analog of the strain parameter 
77 enters 23 is an interesting possibility. 
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VI. Normal stresses and the one- component limit 


The predictions of these theories for homogeneous shear flow in the limit of very large 
77 are compared in Table III. It has been shown 24 that RDT predicts the initial evolution 
of these flows very well. The tensor bij is defined as usual by bij = —Tij/2K + 1/3 6ij. 
The line NLEV refers to the nonlinear eddy viscosity model, Eq. (52) with constant coef- 
ficients C V )C T i,C T 2 -,C T 3 . The unphysical results are shown to motivate the introduction 
of 77 -dependent coefficients; all derivations of this model 4,5 have computed perturbatively 
assuming that the strain rate is small. The results for the explicit algebraic LRR model 
of Pope are shown using the notation of Ref. 25; f >2 and 63 depend on the choice of model 
constants. This model predicts a qualitatively correct shear stress ratio 612 , however the 
normal stress ratios are model dependent. By construction, Eq. (53) recovers the RDT 
normal stresses, the one component limit . 24 It is more important that the present time 
dependent theories, Eqs. (45) and (50) also predict this limit at short times. This fact 
shows once more the importance of the short time corrections to RDT: even if the LRR 
model could be calibrated for agreement with the one component state in the limit of 
large 77 , eddy damping would quickly drive the solution away from this state; the strong 
suppression of eddy damping at short times in the present theory will maintain the RDT 
limit longer. 

The difference between the present theory and LRR can be understood by comparing 
Eqs. (3) and (31): although the production term containing Ai m is treated identically, the 
"rapid pressure strain” term kikpA pm k ~ 2 is treated exactly in the present theory by RDT 
but is modeled in LRR by the terms bilinear in stress and strain rate. Renormalization 
group analysis 27 showed that the LRR model is a rational lowest order approximation for 
the rapid pressure strain term; however, it is not valid for very large strain rates. 

VII. Conclusions 

The present theory should be compared to the LRR stress transport model and to 
RDT. Both LRR and RDT arise from a special choice of the Green’s function, correspond- 
ing respectively to the long time and short time limits of the present theory. The present 
theory reduces to RDT at short times and at large strain rates, and to an eddy damping 
model qualitatively similar to LRR at long times. 
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TABLE I 

Shear Stress Ratios in Homogeneous Shear Flow 
and in Equilibrium Simple Shear Flows 


-r/K 



77 ~ 3.0 

77 ~ 4.4 


Experiments 

.33 

.33 

&& 

Theory 

.32 

.32 


AlgebraicLRR 

.37 

.38 

&&& 

Theory with 




transport corrections 

.32 

.23 

& 


TABLE II 

Stress Ratios in Near-Wall Turbulence 

V -t/K(DNS) 26 -t/K( Theory) AlgebraicLRR —r /if (Corrected Theory) 


6.6 

.178 

.274 


.380 


.246 

8.4 

.159 

.237 


.361 


.195 

9.6 

.149 

.216 


.347 


.170 

10.9 

.138 

.183 


.332 


.139 

12.5 

.127 

.161 


.314 


.117 

14.0 

.116 

.146 


.295 


.105 




TABLE III 




Theoretical predictions 

in the large 77 

limit 




&12 

bn 

b 22 

b 33 



NLEV 

o( V 2 ) 

0(V 3 ) 

0 (,>) 

0(7 7 3 ) 



LRR 

o( v- 1 ) 

4 363 — 62 

4 3&3 4’&2 

8 36 2 



15 36|-fe| 

15 3&1— 

45 36|-6| 



Eq.(53) 

0 ( 77 - 1 ) 

2/3 

-1/3 

-1/3 



RDT 

oiv- 1 ) 

2/3 

-1/3 

-1/3 
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